Models of impurities in valence bond spin chains and ladders 
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We present the class of models of a nonmagnetic impurity in S = i generalized ladder with 
an AKLT-type valence bond ground state, and of a S = = impurity in the S = 1 AKLT chain. 
The ground state in presence of impurity can be found exactly. Recently studied phenomenon of 
local enhancement of antiferromagnetic correlations around the impurity is absent for this family of 
models. 
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Over the last decade, low-dimensional spin systems, 
particularly the Heisenberg spin chains and ladders, 
have contittufid to attract considerable attention of 
researchers. Em The interest to spin ladders is particu- 
larly stimulated by the hope to get some insight into 
the physics of metal-oxide superconductors; in support 
of this hope, superconductivity in the ladder compound 
Sro.4Ca13.6Cu24O41.84 in presence of hole doping and 
high pressure was recently reported^ It is now well es- 
tablished that "regular" (i.e., with only "leg" and "rung" 
exchange couplings) S — h isotropic spin ladders have a 
disordered gapped ground state when the number of legs 
is even, while odd-legged ladders have quasi-long-range 
ordered gapless ground state. On the other hand, "gen- 
eralized" ladders including other couplings can serve as 
interesting toy models with a rich behavior which is often 
very different from that of "regular" models su 

Recently, interesting experimental results on ladders 
doped with nonmagnetic, impurities (Cu substituted by 
Zn) have been obtained:E9 surprisingly, the antiferromag- 
netic (AF) order was found to be stabilized by the doping; 
a similar, behavior has also been observed. in , s pin-Peierls 
chains.til A number of numerical studiesliaO indicated 
that local AF correlations near a nonmagnetic impurity 
are enhanced comparing to the system without vacancies. 
It has been suggested that this phenomenon, as well as 
several other similar effects in one- and two-dimensional 
antiferromagnetsji-j can be explained on a common ba- 
sis using the soi-ealled "pruned" resonating valence bond 
(RVB) pictured Nonmagnetic impurity affects forma- 
tion of instant singlet bonds for spins which are located in 
its immediate vicinity, making some of the bonds geomet- 
rically impossible and thus enhancing the other bonds. 
This explanation is supposed to be rather general and 
does not depend much on the interaction details. 

In this paper I show that for certain models of nonmag- 
netic impurities in generalized S = \ spin ladders with 
exact matrix-product ground states of the type consid- 
ered by us recently,IM local AF correlations are partly or 
completely insensitive to the presence of impurity. 

Consider the model of a vacancy in the generalized 
S = 5 ladder with additional diagonal and biquadratic 
interactions, described by the following Hamiltonian: 



H — ^ hi,i+i 



(1) 



Si,i+rS2,i+i) 



+ j£ (Sij-Sij + S2i'S2j) + Jrj(Sli-S2j + Sgj'Slj) 
+ V&tSii-SyXSai-Siy) + l ; '/ ; ,(S;,.S 2 ,|(S,,.S;,: , 

here the indices 1 and 2 distinguish lower and upper legs, 
and i labels rungs (see Fig. |lj), and the terms involving 
the vacancy site S*2,o are implicitly assumed to be missing 
in /io,i and h- 



J 



-1,0- 

= 1, and V, 



The "bulk" couplings Jr, J l £ l+1 



= v, 



D — a., auu v LL - v DD = I do not depend on i, 
Jr is a free parameter, and we have introduced the extra 
"edge" interaction between the rungs —1 and 1 across 
the vacancy to make the problem solvable. 

In absence of the vacancy the model ([!]) describes the 
generalized Bose-Gayen model as introduced in Ref. ||, at 
the special value of the leg/diagonal coupling ratio equal 
to 1. At Jr > I its ground state is a product of singlet 
bonds along the ladder rungs, and we will be interested 
in the interval Jr < §, where the ground state coincides 
with that of the effective S = 1 Affleck-Kennedy-Lieb- 
Tasaki (AKLT) chain,B whose S =1 spin^are formed by 
the triplet degrees of freedom of the rungs£i This effective 
AKLT ground state can be conveniently written in a form 
of the so-called matrix product state:E30 



*o = tr(Q ffl ), 9i=-j= 



\to)i 
V2\t. 



-y/2\t+) 
-|*o> 



(2) 



where |£ M )i, = 0, ±1 are the triplet states of the i-th 
rung. The ground state energy per rung isa 

E = -13/10 + J R /A. 

We will look for the wave function of the ground state 
in presence of the impurity in the form of the following 
matrix product: 



*o" P = tT (g-N ■ ■ ■ g-iGogi ■ ■ ■ ffjv) , 



(3) 



where the matrix Go corresponding to the unpaired spin 
at the 0-th rung is chosen from the requirement that 'J™ 13 



1 



has both the total spin and its_j2-projection equal to |; 
the most general form of Go isEa 
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x being a free parameter. Physically, the wave function 
'J™ 13 describes a superposition [x / y/S)^^, 



1,1/2 



,0,1/2 
1/2 ^ *l/2 

where \H lad ' 1 / 2 denotes a wave function with the total 

Jtot 

spin jtot composed from the states of the unpaired spin 
\ and the states of the rest of the ladder having to- 
tal spin jiad! in ^1/2 the unpaired spin is completely 
decoupled from the rest of the ladder (which is in the 
effective AKLT state with one valence bond across the 
impurity), while in ^!vV 2 it is coupled with the edge 

Kennedy- Tasaki tripletEj into the state with j tot = i, 
see Fig. |[ 

Further, the Hamiltonian (|l]) conserves parity with re- 
spect to the mirror transformation i t— > —i, and one can 

see that ^jy 1 / 2 and ^J/ 1 / 2 have different parities. The 
solution with completely decoupled unpaired spin is not 
very interesting, so that we will look for the ground state 
wave function of the form (||), (0) with x = 0. Following 
the approach outlined in Ref. |7l we demand that the 
local Hamiltonian /ii mp , defined as (recall that the terms 
with 52,o should be dropped) 



h-lfi + ho i + /l-H — £q : 



(5) 



where eg is a free parameter, annihilates all states con- 
tained in the matrix product <?_iGo<?+i, and that all 
other eigenstates of h- lmp have positive energies. These 
conditions are sufficient for ^ mp to be the ground state 
of H. The- construction routine is well described in 
literature JlJ'Q'EI so that I only briefly address it here. 

The states of the [—1,0, 1] block can be classified into 
multiplets $fj m , where j is the total spin of the block and 
to is its z-projection. In total, there are ten multiplets 
(five with j = \, four with j = |, and one with j = |); 
one can however straightforwardly check that the matrix 
product g_iGo5+i contains only states of the following 
three multiplets: 



i,m> 



(6) 



here ^ Sb ' Sab denotes the state of the [-1,0, 1] block 
with the total spin j, Sa, Sb, and Sab being the total 
momenta of the —1th, -l-lth rung and the [—1, 1] block, 
respectively. The local Hamiltonian h lmv should annihi- 
late the states (@), so that it can be generally written 
as a projector onto the subspace of the remaining seven 
multiplets, 



^ =^(V 1 i 01 



*? 2 =t/j 

— m I 



000 
. m '■ 



(7) 



ttS 1 = 






*3 2 = 




^ e 3 ' 3 = 






*5 = 

2' m 


iPl 12 



We make a further simplification, assuming that hi mp 
does not mix the above multiplets, so that 



E E E A 



(8) 



where all A"-^ should be positive to ensure that (|) is the 
ground state. Demanding further that this structure is 
compatible with the particular form of the Hamiltonian 
([j]), one arrives at the following family of solutions for 
the coupling constants and the parameter Sq: 



J 
J 



L 1J = x[)\/2+(l + J R )/4, V-t 1 
_ n _ \(i) 



Jr. + 2A^V 2 — 1 



-2A 



(i) 

1/2 



-19/16+ Jr/A, 



(9) 



where A^V 2 plays the role of a free parameter, and the 
expressions for the other eigenvalues are 



Xy 2 = 1 — Jr, A^2 — 1/4 - (A^2 + Jr)/2, 
A« =3/2 + Af; 2I Af; 2 =3/2, (10) 
A ( 3 ) -9/«:_ AW 



>(3) 



(1) 



v 3/2 



2/5 - (Ai% + Jr)/5, A 5/2 = 5/2. 



The parameter eo has the meaning of a ground state en- 
ergy of the [—1, 0, 1] block, thus the states with and with- 
out a vacancy differ in energy by the value eo — 2£'o- The 
conditions of positivity of h- lmp require that 



J ^2 



A 



(i) 

1/2' 



(11) 



The most symmetric solution from the above family is 
achieved by setting A^ 



(i) 

/2 



(1 - 2J R )/A, J R < 1/2, then 

J- M = J" 1 ' 1 = 3/8 and V^ X = V D ^ = -1/2. _ 

Using the standard matrix product technique, it is easy 
to calculate the spin correlation functions and distribu- 
tion of the excess spin in the state (||) . The mean value 
of S z at each site is given by 



Kl-i 



(12) 



here |«| > 1. Following Ref. [1^, we calculate the spin cor- 
relation functions along the ladder legs, with the starting 
site being next to the vacancy, and compare them to the 
correlations in absence of the vacancy. Quite surprisingly, 
one finds that the AF correlations are not at all affected 
by the presence of a vacancy: 
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(13) 



here a = x, y, z, "w/o" in the second line means "without 
vacancy," and |i| > 1. Note that despite the presence of 
the excess spin, the spin correlations remain isotropic. 

One can show that the above features (insensitivity 
of AF correlations to the presence of a vacancy and 
their isotropic character) survive also in more compli- 
cated matrix-product-solvable models: the ansatz (^), 
(||) can be obviously used in its most general form, with 
the "bulk" matrices gt iwcJuding singlet degrees of free- 
dom of the ladder rungs ,E3 



9i i-> gi(u) 



I 



V3 



u\s)i + \to)i 
V2\t-)i 



-V2\t+)i 
u\s)i - \to)i 



then generally the Hamiltonian (|l|) will not be invariant 
under the parity transformation i — > — i, so that the ma- 
trix (Q) can be also used in its general form with arbitrary 
x. [The corresponding family of impurity models with ex- 
act ground states can be obtained for this general ansatz, 
exactly in the same way as we did above; however, the 
resulting model Hamiltonians are extremely cumbersome 
and therefore we do not present those solutions here.] For 
such a state, the above formula for the correlations will 
change as follows: 



Wl,l°M 



\°l,0°l,i 



\°2,1 °2,i 



(^2,n^2,n+i)w/o — 1 % I (u 2 + 3) . 

u 2 -l 



(14) 



(l + x)q 



u 2 + 3 



One can see that the AF correlations along the legs are 
not affected, except for the correlations involving the un- 
paired spin Si^o, which are enhanced for u being in the 
interval between — x and [x — 3)/(l + x) and suppressed 
otherwise. In valence-bond-type models the decay of all 
correlations is purely exponential for all distances, and 
presence of the impurity can only change the prefactor in 
front of the exponent; accidentally, for the chosen ansatz 
(||,f|) the changes coming from the excess spin and "dis- 
tortions" due to the presence of a vacancy completely 
compensate each other. The spin excess distribution is 
also modified and is generally asymmetric: 



2(^-1)^1 2(x-v)qW 



(1 - <7U)(3 + X 2 ) 
(x 2 - 1) 



(u + a){3 + x 2 Y 



2(3 + x 2 )' 



cr = sgn(«), |«|>1. 



Finally, one can observe that the model of a vacancy 
in the S = ^ ladder can be reformulated as a model of 
the S = \ impurity in the 5 = 1 AKLT chain. Consider 
the model described by the following Hamiltonian: 



H - }Xhf^f + h^f^t) + /limp, 



(15) 



%>i 



h imp = (J+Si + J_S_i) • t + J'(S_i • Si) - e 
+(S_i • Si){(^ + Si + F_S_i) • r} + V"(S_i • Si) 2 , 

here h^f LT = Si ■ Sj + 1 (S t • S 3 ) 2 is the local Hamiltonian 

of the AKLT chain in the bulk, and h lmp describes the 
interaction induced by presence of the impurity spin r, 
see Fig. f| the parameter £o being just a constant energy 
shift having the meaning of the ground state energy of the 
[— 1, r, 1] block. Using the ansatz (||), (Q), one can repeat 
the entire construction routine as described above for the 
ladder, and obtain the following family of Hamiltonians 
for which "J™ 13 is the exact ground state: 

J± = {5(3 - x 2 )/9 ± 2x}\ 3/2 + 5A 5/2 /9, 

= -(5 + ^ 2 )A 3/2 /3 + 5A 5/2 /6, (16) 
V± = {-5(3 + x 2 )/9 ± 4x/3}A 3/2 + 5A 5/2 /9, 
^ = -(15 + a; 2 )A3 /2 /9 + 5A 5/2 /18, 
£ = (30-2x 2 )A 3/2 /9 + 5A 5/2 /9. 

Here A 3 / 2 > 0, A 5 / 2 > are the eigenvalues of /ii mp 
corresponding to the multiplets 



y. For x ^ 0, 00 the Eqs. (pf) 



! = V" 2 , 

' 2' m 



respectively. For x 7^ 0, 00 the Eqs. ( pL6[ ) describe mod- 
els with an asymmetric impurity. Again, as in case 
of the ladder, one can straightforwardly check that the 
"edge" spin correlation function in presence of the impu- 
rity (S±S") just coincides with that in the bulk, indepen- 
dently of the value of x. The above family contains two 
interesting solutions: one is achieved by setting x = 0, 
A5/2 = 3A 3 / 2 and describes the simple symmetric model 
without biquadratic terms involving the impurity spin r: 



J± = J, J' 



-V 



1 



:J, V±=0. 



(17) 



Another solution corresponds to A 3 / 2 = 0, then the 
ground state of the model is twofold degenerate since 
both even and odd-parity wave functions ^^(x = 0) 
and ^™ lp (x = 00) are eigenstates with the same energy. 
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FIG. 1. Nonmagnetic impurity in the generalized S = \ 
spin ladder as described by the Hamiltonian (jl|), V's denote 
the biquadratic couplings. 
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FIG. 2. Schematic valence-bond representation of the 



wave functions ^l'^ 2 an( A ^1/2^ contained in (|3|); solid and 
dashed lines denote singlet and triplet valence bond links, re- 
spectively. Solid ovals indicate that spins on each rung are 
coupled into effective triplet; dashed oval in the bottom pic- 
ture denotes that the triplet valence bond and the unpaired 
spin are coupled into a spin-^ state. 
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FIG. 3. S = \ impurity in the S = 1 AKLT chain as 
described by the Hamiltonian (|l5|), r is the impurity spin, 
and V^'s indicate the biquadratic couplings. 
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